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our previously published results for the nonadiabatic contributions we show that the adiabatic 
contribution, which is proportional to T 2 at low temperatures and goes as T 3 at high temperatures, 
dominates the nonadiabatic contribution for temperatures above a cross-over temperature, T c , which 
is between 0.5 and 0.8 T m , where T m is the melting temperature of the metal. The nonadiabatic 
contribution falls as T _1 for temperatures roughly above the average phonon frequency. 
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I. INTRODUCTION 

The crystal free energy consists of contributions from 
the static lattice potential, the phonons, the electronic 
excitations, phonon-phonon interactions, and electron- 
phonon interactions (Wallace^, eq. (20.1)). The last 
term is generally the smallest, but little is known about 
the size of this contributions except at very low tem- 
peratures. We therefore study the electron-phonon free 
energy, F ep , in order to assess its contribution to the 
thermodynamic properties of metals at all temperatures 
to melting. 

The total electron-phonon free energy was derived by 
Eliashberg— in a paper on superconductivity. In this for- 
mulation, the electron-phonon interaction in the crys- 
tal ground state is double counted. The ground state 
subtraction, which corrects for the ground state double 
counting, was derived in Wallace^, eq. (25.26). It is well 
known that electron-phonon interactions contribute to 
the low temperature electronic specific heat C e i , accord- 
ing to 



C el =T bs (l + A) T. 



(1) 



The bare electron contribution is r^T, where is pro- 
portional to the electronic density of states at the Fermi 
level determined by bandstructure. The nonadiabatic 
electron-phonon contribution gives the term in A, and 
can be quite large, e.g. A as 1.5 for lead^. In the early 
years, many theoretical calculations of A were carried out 
for the simple metals, showing rather good agreement 
between theory and experiment. These results were re- 
viewed by Grimvall 6 , Tables III- VI; see also Wallace 3 , 
Tables 27 and 28. Grimvall 6 evaluated the Eliashberg 



formula for lead, and concluded that F ep vanishes at tem- 
peratures above the phonon characteristic temperature. 
This conclusion will be revised in the present report. 

The adiabatic approximation^ rests on the expansion 
of the coupled nuclear-electron Hamiltonian in powers of 
m/M , the ratio of electron to nuclear mass. When terms 
of order m/M are neglected, the electronic wavefunctions 
do not see the nuclear motion, and depend only on the 
static nuclear positions^. The case of lattice dynamics is 
treated by Born and Huang 8 , p. 171. In the partition 
function, the same property is revealed by treating the 
operation of the nuclear kinetic energy on the electronic 
wavefunctions as a small effect 9 . When this operation 
is neglected entirely, which is appropriate at high tem- 
peratures where the nuclear motion is classical, only the 
adiabatic part F ad survives in F ep (Wallace^, p. 91- 
93). It is therefore seen that the electronic states are 
not mixed in the adiabatic part, and also that F ad is the 
dominant part at high temperatures. The nonadiabatic 
part Fi a + Fg 0, arises from mixing of electronic states by 
the nuclear kinetic energy, and becomes unimportant at 
high temperatures. 

Allen and Heine^S separated electron-phonon effects 
into adiabatic and nonadiabatic contributions, and stud- 
ied the adiabatic part of the electron energy shifts. For 
the Eliashberg formulation, Allen and Hui 11 studied the 
adiabatic contribution to the high temperature specific 
heat. Present results will be compared with Allen and 
Hui in Sec. HTH The leading correction to the linear tem- 
perature dependence in the specific heat given in eq. 
is a nonadiabatic term which goes as T 3 lnT. This has 
been studied extensively in the literature starting with 
Buckingham— and Buckingham and Schafroth— . See 
reference^ ancU^ and references therein for more recent 
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work. 

We have previously addressed the problem of calculat- 
ing the contribution to the free energy from the electron- 
phonon interaction 5 -. We found the nonadiabatic part of 
this contribution, F™, to second order in the interaction. 
We calculated F% a for four nearly-free-electron metals, 
Na, K, Al, and Pb, for temperatures between zero and 
roughly 1.5 times the melting temperature. 

As in the calculation of the nonadiabatic contribu- 
tion, our calculations are done for a constant density 
p to eliminate concern for the density dependence of 
phonon frequencies and elect ron-phonon interaction ma- 
trix elements. The density is that at the temperature T ( „ 
where the phonon frequencies are measured. The melting 
temperature at this density is higher than the custom- 
ary zero-pressure melting temperature. Our calculations 
cover the range from T — to above T m . 

The paper is organized as follows. In Sec. [H] we 
present the expression for the adiabatic and nonadiabatic 



parts of F ep , and we discuss the formulation for metals 
in general, the simplification for the nearly-free-electron 
metals studied here, and the groundstate subtraction. 
Asymptotic temperature dependences are derived for the 
adiabatic part. In Sec. IHII our results are presented with 
a detailed discussion of the numerical methods used, the 
adiabatic and nonadiabatic parts are compared, and our 
results are compared with previous work. Our conclu- 
sions are summarized in IIVI 



II. ADIABATIC CONTRIBUTION TO THE 
ELECTRON-PHONON FREE ENERGY 

A. Analytic Form of the Free Energy 

The electron-phonon contribution to the free energy 
can be written in three pieces, F ep = F ad + F™ a + F 2 na , 
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Our results are calculated and quoted per atom. 
Throughout this paper we will use the following nomen- 
clature: fp is the Fermi-Dirac distribution function at fi- 
nite temperature and gp is the same at T = 0. is the 
Bose-Einstein distribution function at finite temperature 
and is the polarization vector of the phonon branch 

A for wave vector k which is inside the Brillouin zone. Q 
is a reciprocal lattice vector and loj: x is the frequency of 
a phonon mode. U(k + Q) is the Fourier transform of the 
pseudopotential for momentum transfer k + Q 

To remind the reader of the physical meaning of the 
three terms in eqs. 10 - HJ ( a more detailed discussion 
can be found in Bock et al. 5 ) a quick summary: Eq. J2J, 
F ad , expresses the thermally averaged vibrational contri- 
butions to the excited electronic energies. Eqs. Q and 
Q, F™, describe the non-adiabatic corrections to all 
electronic energy levels and take into account the mixing 
of electron states due to the ion motion. 

In the free energy formulation, the electron-phonon 
interaction is treated in second order perturbation 



theory 1 ! 2 ! 3 ! 4 ! 5 ! 6 ! 10 ! 11 ! 14 ! 15 . For a general metal, the elec- 
trons are presumed to have band structure, and the en- 
ergy denominators are band electron energies Eg. For 
nearly-free-electron metals, band structure effects may 
be treated in pseudopotential perturbation theory, where 
in zeroth order the electron energies are the free electron 
energies eg. This is why free electron energies appear 
in the denominator of eqs. J2J~0- In these equations, 
the band structure effects, to second order in the pseu- 
dopotential, are contained in the last term in brackets in 
eq. J2J). Pseudopotential perturbation theory has been 
extensively developed over many years, and pseudopo- 
tential parameters have been calibrated to experimental 
data such as equilibrium density and bulk modulus. We 
use these calibrated pseudopotentials here, so that our 
models have no free parameters. 

Finally we need to clarify the problem of groundstate 
double counting. The potential for the nuclear mo- 
tion, the "adiabatic potential" , is precisely the electronic 
groundstate energy as a function of static nuclear posi- 
tions. Since this potential has been put into the phonon 
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Hamiltonian, the electronic groundstate energy has to be 
subtracted from the electron-phonon Hamiltonian. In 
this way the electronic statistical mechanics expresses 
only electronic excitations from the groundstate. The 
groundstate subtraction is an adiabatic effect, and is ex- 
pressed by the term (—gp) in eq. (0. The decomposition 
of the total Hamiltonian for a metal crystal, and deriva- 
tion of the corresponding free energy, including eqs. 
is given in Wallace^. 

B. Analytic temperature dependence 

The temperature dependence of F ad , eq. (J2J, arises 
from the product of the Fermi-Dirac factor (fp — gp) and 
the phonon factor (n^ A + 1/2). The value at T = 0, i.e. 
the constant term in F ad , vanishes because (fp — gp) van- 
ishes at T = 0. Using a Sommerfeld expansion (e.g. pp. 
45 in Ashcroft and Mermin— ) of eq. (J2J) , the Fermi-Dirac 
factor gives a quadratic temperature dependence in lead- 
ing order, plus higher order terms which are of relative 
order (kgT/ep) 2 and so can be neglected at temperatures 
to well above T m . The phonon factor reduces to the con- 
stant 1/2 at very low temperatures, i.e. at fcgT <C the 
mean phonon energy, (tjio), so that the net dependence 
at very low temperatures is 

F ad = B 2 {k B Tf . (5) 

At higher temperatures, i.e. at UbT > (Hu>), each phonon 
factor can be expanded as 

{^ + 2)-^ + ^ + -^ + - (6) 

Hence the temperature dependence at k B T > (hui) be- 
comes 

F ad = B 3 (k B Tf + ■■■ (7) 

As we will show in a later section, our numerical results 
confirm this analytic temperature dependence. 

III. RESULTS 

A. Phonon Model 

When we calculated the non-adiabatic contributions 
we found that we could replace the phonon dispersion 
with a constant frequency with high accuracy which sim- 
plified our numerical evaluation. In fact, we found that 
the result of our calculation was not sensitive to the value 
of this frequency, u>e, and that an Einstein model gives 
the correct temperature dependence. We reasoned that 
this insensitivity is due to the fact that the integrands 
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FIG. 1: (Color online) ^-dependence of F ad for Na at T = 
200 K. 

of F™2 are finite and well-behaved as a function of wg v 
The adiabatic contribution on the other hand diverges as 
lue — > as illustrated in Fig. ^and is therefore very sen- 
sitive to the choice of lie- Although the Einstein model 
will give the correct temperature dependence at low and 
high temperatures as listed in eqs. |JSJ and we have no 
way to find the Einstein frequency which will give us the 
correct magnitude without doing a calculation with the 
full phonon dispersion. This is a serious problem as can 
be seen from the extreme dependence on loe as shown 
in Fig. ^ Since in the case of F ad the calculation using 
the full phonon dispersion is not prohibitive, we chose for 
this report to use a fully fc-dependent phonon dispersion. 

In the following we calculated the phonon dispersion 
and eigenvectors from Born-von Karman force constants 
to get a more realistic representation of the real phonon 
spectrum with its three branches (the force constants 
were taken from Dederichs et al.— £). 

B. Numerical Techniques 

Eq. is written in terms of sums over the electron 
and phonon momenta. We are using a free electron dis- 
persion spectrum which is isotropic in the electron mo- 
mentum. This combined with the fact that |p| is bounded 
by (fp — gp) from above and below to |p| « p F to within 
a few ksT makes it numerically more convenient for us 
to rewrite J2p as an integral, 




The sum over phonon momenta however is evaluated 
more conveniently as a sum because of the two terms 
in the curly brackets in F ad . The momentum transfer 
k + Q enters the second term only with the reciprocal 
lattice vector Q. Were we to integrate over the final elec- 
tron momentum, p' = p + k + Q, the second term would 



introduce discrete steps into the integrand as we change ming over k and Q instead converges much faster. The 
p'. This is very difficult to handle numerically. Sum- expression we evaluated is given by 
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(9) 




where V denotes the Cauchy principal value. 

At this point we would like to address the question of 
convergence in our calculation since the sum on k + Q 
is not obviously bounded by an upper limit. We had 
found for the two nonadiabatic parts strong convergence 
as a function of p' due to a combination of pseudopoten- 
tial factor and the energy denominator. In the case of 
the adiabatic contribution however, this is not quite so 
obvious. The energy denominator is only linear in the en- 
ergy difference and we therefore expect the convergence 
of the two single terms in F ad to go as Q~ 2 , which is 
quite weak. We do suspect however that the two terms 
in the curly brackets will exhibit cancellation to some 
degree and might improve the convergence behavior of 
F ad . In Fig. |21 we plot the two terms separately and 
combined as a function of magnitude of Q. Despite the 
slow convergence behavior of the terms by themselves, 
their combination converges strongly. 



C. Comparing the adiabatic with the nonadiabatic 
results 

As mentioned earlier, we calculated F ad using force 
constant models for the phonon dispersion for the ele- 
mental metals Na, K, Al, and Pb. Our result for Na is 
shown in Fig. [3| A polynomial fit of the temperature de- 



pendence in the low- and high-temperature regime shows 
that our initial analysis was correct and F ad oc T 2 for low 
temperatures and that it turns into oc T 3 for higher tem- 
peratures. 

The sign of F ad for all metals studied here except Al 
is positive. It is not obvious from eq. © why this is the 
case. We can speculate that the sign is determined by 
the relative magnitudes of the two terms in parentheses 
in F ad since we know from our convergence study (Fig. 
EJ) that the two terms are roughly of equal strength but 
of opposite sign. The difference between the two terms is 
the scattering momentum which is the full phonon mo- 
mentum in the first term (k+Q) and the reciprocal lattice 
vector only in the second term (Q). 

The Fermi-Dirac factor (fp — gp) ensures that the first 
electron momentum, p, is essentially confined to the 
Fermi surface, pf- For the denominators to be small, 
the second electron momentum, (p + k + Q) and (p + Q), 
has to therefore lie close to the Fermi surface. Although 
the exact shape of the Fermi surface can be quite compli- 
cated we would like to restrict this argument to a spher- 
ical Fermi surface since we are using the pseudopotential 
formalism in this calculation and do not take separate 
bands into account. In this case it is plausible to spec- 
ulate that a change of the radius of the Fermi surface 
might be able to effect an overall sign change in F ad . 

We explore this possibility by means of a "simplified 
F ad " , given by eq. © with the pseudopotential set to 



5 



200 
150 
E 100 
2 50 

a o 
I - 50 

I -100 
-150 
-200 

















Na — i — 
Al 






n /Ah - 


1 .3436 A" 1 








Pi (Na) 


= 1.0503 


V 1 


P2 (f 


Ja) = 

r. 


1.4854 A" ■ 

Al\ 1 CCiA A"' 


- 

- 










J 

i 


f 




™/ - 1 ■ 

< 






> 




Mai 


= 0. 


►2 A " 1 




+ 






IV 


= 1:75 A" 1 




































.8 1 


1.2 


.4 


1.6 


1.8 





ficticious p F [A" 



FIG. 4: (Color online) Dependence of sign of F ad as a func- 
tion of p F for Na at T = 200 K and for Al at T = 500 K. For 

illustration purposes, the lengths pi,2 = 1/2 |Qi,aj of the first 

two reciprocal lattice vectors are shown. 



a constant and with a "fictitious pf" that is allowed to 
vary. This amounts to pretending that we can change 
the number of free electrons, hence change the size of 
the Fermi surface, without changing anything else. In 
this case the essential part of the second term in eq. © 
is the p-angle average of (e^— e^.^) -1 with \p\ equal to 
the fictitious pf and with Q any reciprocal lattice vector. 
This integral diverges as the fictitious pf passes through 
p n = 1/2 Q n , where Q\ is the magnitude of the first 
reciprocal lattice vector, and so on for Q 2 , etc. At the 
same time, the first term in eq. is well behaved as 
a function of the fictitious pp- Numerical evaluation of 
the simplified F ad is shown in Fig. 0] for Na and Al, 
and reveals the expected logarithmic discontinuities at 
p n = 1/2 Q n . Fig 0] supports the view that the sign of 
pad d e p en d s directly on the size of the Fermi surface 
relative to the Brillouin zone surface. 

Our results for Na for F ad including our previous re- 



sults for F™ are shown in Fig. [SJ At low temperatures, 
the adiabatic contribution vanishes and the nonadiabatic 
contribution approaches a constant. In the low tempera- 
ture regime, the free energy is dominated by the nonadi- 
abatic contribution. In the high temperature regime, the 
adiabatic contribution increases as T 3 and the nonadia- 
batic contribution slowly vanishes. The adiabatic contri- 
bution dominates in this temperature regime. At about 
sa 320 K there is a cross-over between the two contribu- 
tions. For reference we included the temperature of an 
average phonon frequency and the melting temperature 
of Na in the graph. The cross-over temperature is well 
between either one of these two temperatures. 

Fig. shows our results for K. We find the same qual- 
itative behavior. The cross-over temperature is always 
well between the two reference temperatures. The same 
is seen for our results for Pb, shown in Fig. [7| 

In the case of Al we repeated the calculation for both 
pseudopotential models. Our results are shown in Fig. [S] 
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FIG. 5: (Color online) F ad and F 2 " a for Na. The cross-over 
temperature is shown. 
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FIG. 6: (Color online) F ad and F™ for K. The cross-over 
temperature is shown. 



The different pseudopotentials effect a slightly different 
curvature of the temperature dependence in the adiabatic 
contribution. As we found previously, the nonadiabatic 
contribution is shifted in energy by about 5 x 10~ 5 eV. 
The net result of these two effects on the cross-over tem- 
perature (between the magnitudes of F ad and F™) is 
to shift this temperature by only roughly 20 K. The dif- 
ference between the two pseudopotential models is pre- 
sumable well within the accuracy of the pseudopotential 
method and our calculation itself. It does not make any 
difference which model is chosen. 

All of the cross-over temperatures are summarized in 
table □ 

In order to assess the overall importance of the adia- 
batic electron-phonon free energy at high temperatures, 
we compare it with the bare electron free energy F e i . The 
leading order Sommerfeld expansion of F e i , which is quite 
accurate to T m for the nearly-free-electron metals, is 
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TABLE I: Cross-Over temperatures and other quantities 
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FIG. 7: (Color online) F ad and F? a for Pb The cross-over 
temperature is shown. 
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FIG. 8: (Color online) F ad and F? a for Al. Note that F ad 
is negative but drawn here as —F ad . The cross-over temper- 
ature is shown. Both pseudopotentials are used. The inset 
shows the cross-over region. The datasets are labeled as in 
the larger plot. 



Since F ad increases approximately as T 3 at high T, the 
ratio F ad /F e i has its maximum magnitude at T m . This 
ratio is listed in Tabled based on our calculations and the 
electronic density of states of Moruzzi et al.™. Our cor- 
responding values for S ad /S e i at T m , and C ad /C e i at T m , 
are also listed in Tabled These results allow us to make 
estimates for nearly-free-electron metals, of whether or 
not the adiabatic electron-phonon contribution is signif- 
icant in a given property of a given metal. 



D. Comparison with previous work 

Previous work on the electron-phonon contribution 
was done on the entropy instead of the free energy. The 
thermodynamic relation 
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FIG. 9: (Color online) The free energy contributions F and 
F 2 na for Pb. 



relates the two quantities. In order to see how our re- 
sult for the free energy affects the entropy, we took the 
derivative of the adiabatic and nonadiabatic contribu- 
tions and plotted the entropy for temperatures between 
100 K and 1000 K for lead in Fig. ED The nonadia- 
batic entropy slowly vanishes with increasing tempera- 
ture, whereas the adiabatic contribution increases. The 
electron-phonon contribution to the entropy in lead was 
calculated by GrimvalU^. We found previously that for 
temperatures up to T < 1.4 Te the nonadiabatic con- 
tribution to the free energy is sufficient to achieve good 
agreement with GrimvaH's entropy calculation. We see 
now why. The adiabatic contribution is too small to con- 
tribute noticeably to the total entropy. At higher tem- 
peratures however, this will not be the case anymore and 
the adiabatic contribution will become important. As we 
see from Fig. E3 S increases with temperature and at 
around T — 240 K \S ad \ > IS^I- Grimvall's conclusion 
that S ep (T — > oo ) = is therefore not correct. We also 
notice that, maybe not surprisingly, the cross-over tem- 
perature for the free energy is higher than the cross-over 
temperature for the entropy. 

Allen and Hui 11 studied the adiabatic electron-phonon 
specific heat as derived from the Eliashberg free energy. 
They argued that for A15 metals, which have narrow 
peaks in the electronic density of states, the adiabatic 
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FIG. 10: (Color online) The entropy contributions S ad and 
for Pb. The free energy is shown for comparison in Fig 

m 

contribution can be as large as the bare electronic spe- 
cific heat at high temperatures. Our results in Table 
H] are not inconsistent with this possibility. Differences 
in our procedure and theirs prevent a further compar- 
ison of results. First, their electron-phonon free en- 
ergy does not contain the groundstate subtraction. Wc 
note, however, that while the groundstate subtraction 
contributes to the free energy and entropy at all temper- 
atures, its contribution to the specific heat vanishes at 
high temperatures. Second, our adiabatic specific heat is 
C ad _ -T (d 2 F ad /dT 2 ) v , while Allen and Hui remove 

certain terms from this C ad , to be placed in the bare 
phonon and bare electron specific heats. This, of course, 
is not an error but a matter of choice. Nevertheless all 
the terms in question arise from eq. J2J for F ad and are 
therefore genuine electron-phonon interaction terms. 

IV. CONCLUSIONS 

A. Structure of the Free Energy 

The exact free energy of a crystalline elemental metal 
can be expressed asi^ 

F = $o00 + F el (V, T) + F ep (V, T) 

+ F ph (V,T) + F anh (V,T), (12) 

where $o is the static lattice potential, F p h represents 
quasiharmonic phonons, F an h represents phonon-phonon 
interactions, F e i represents static-lattice electronic exci- 
tations, and F ep is the remainder of the free energy ex- 
pressed in terms of interactions between electronic excita- 
tions and phonons. F ep is the smallest and most compli- 
cated term in eq. (|12|) . I n leading order thermodynamic 
perturbation theory, F ep is composed of a nonadiabatic 
part, F na = F[ la +F£ a , and an adiabatic part, F ad . The 
nonadiabatic part arises from the mixing of electronic 



states due to the motion of the ions, and strongly affects 
the lowest electronic excitations. Hence F na dominates 
F ep at low temperatures. The adiabatic part, F ad , ex- 
presses the averaging over the thermal motion of the ions 
of each electronic excited state, and is formally the lead- 
ing contribution to F ep at high temperatures where the 
ion motion is classical. Hence there exists a crossover 
temperature T c , below which the major contribution to 
F ep is F na , and above which it is F ad . 



B. Sign and Magnitude Estimates 

The following properties are characteristic of our re- 
sults for Na, K, Al, and Pb. We suppose these properties 
are common but not without exception among metals in 
general. 

1. F? a is negligible compared to F 2 " a for < T < T m . 

2. F 2 a is positive at all T. From a positive value at 
T = 0, F 2 a decreases as T 2 at first, then decreases 
ever more slowly with increasing T (Figs. 4, 6, 8, 
9 of Bock et alA). 

3. At very low temperatures, F™ is the only signifi- 
cant contribution to F ep , and has the form Kj 10 = 
C 2 + A 2 T 2 , where C 2 > and A 2 < 0. C 2 consti- 
tutes an entirely negligible contribution to the elec- 
tronic groundstate energy. The temperature de- 
pendence is the same as F e i = —(l/2)Tb s T 2 at low 
temperatures, so that F 2 a causes the well-known 
electron-phonon correction to the bare electronic 
specific heat. We have C e i + C ep = (Tb s — 2A 2 ) T, 
where |2A2/Fb s | ranges from a few percent to 
around 2. 

4. F ad can be of either sign. F ad is negligible com- 
pared to F 2 a at low temperatures, but |F ad | > F 2 a 
at T > T c . This is because of the strong temper- 
ature dependence of F ad = B 3 {k B T) 3 at k B T > 
(fku). Though |i^ aci /F| la | increases with tempera- 
ture, F 2 a it not always negligible at T m (Figs. 0- 
[SJ. On the other hand, S 2 a is negligible compared 
to \S ad \ at high temperatures (Fig. HU|) . 

5. The ratio |F ad /F e ;| reaches its maximum at T m , 
where it ranges from a few percent to 0.2 (Table 
IH) . The ratio \S ad /S el \ is larger, and \C ad /C e[ \ is 
larger still (Table QJ, so that the adiabatic contri- 
bution to the specific heat can become important 
at high temperatures. 

6. The theory is a complex mixture of electronic ex- 
citations and phonons, hence the crossover tem- 
perature T c does not scale with either ep or {fna). 
From our calculations, ksT c is well above (fuv), and 
T c /T m is in the range 0.5 - 0.8. 
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7. Based on our results shown in Fig. we expect 
to see a significant dependence of F ad on the con- 
centrations of the constituents in a real system in 
which the number of electrons in the conduction 
band can be controlled, as in an intermetallic com- 
pound for instance. 

It is possible to carry out an accurate numerical calcu- 
lation of the low-temperature nonadiabatic coefficient A2 
(see e.g. eq. (25.83) of Wallace 3 ). A good estimate can 
also be obtained from an Einstein approximation with uje 
determined from Hoje = {hx>) (Bock et al^). As shown 
here (Sec. MI A(l . an accurate calculation of F ad requires 
the use of realistic phonon frequencies and eigenvectors. 



While an Einstein model does give the correct tempera- 
ture dependence of F ad , it does not give a reliable mag- 
nitude (eqs. (28.44) and (28.52) and line 3 of Table 27 
of Wallace 3 ). Based on a numerical evaluation for Pb to 
around 100 K, Grimvall 6 concluded that S ep is negligible 
for temperatures above the mean phonon energy. Our 
results show that this is not the case however, since at 
higher temperatures the adiabatic contribution starts to 
dominate and the total entropy rises, shown in Fig. ^] 
Allen and Hui— argued that C ad can become as large as 
C e i for certain metals at high temperatures. Our results 
for nearly-free-electron metals are not inconsistent with 
such behavior. 
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